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Abstract. These notes describe some applications of the analysis of ordinary differen- 
tial equations to the study of the viscous approximation of conservation laws in one space 
dimension. The exposition mostly focuses on the analysis of invariant manifolds like the 
center manifold and the stable manifold. The last section addresses a more specific issue 
and describes a possible way of extending the notions of center and stable manifold to some 
classes of singular ordinary differential equations arising in the study of the Navier-Stokes 
equation in one space variable. 



1. Introduction 

In recent years, techniques coming from the study of ordinary differential equations (ODEs) 
have been applied to the analysis of systems of conservation laws in one space dimension, 

u t + f(u) x = where u(t,x) E R N and (t,x) G [0, +00 [xR. (1.1) 

In the previous expression, t and x denote the partial derivatives with respect to the variables 
t and x respectively and the flux / : M. N — > JBL N is a function of class C 2 . In view of physical 
considerations (see e.g. Dafermos [7j) it is of great interest considering the second order 
approximation 

ut + f(u% = e(B(u £ )ul) x , (1.2) 

where B is an N x N "viscosity matrix" satisfying suitable conditions and e is a parameter, 
e — > + . In particular, the Navier-Stokes equation in one space variable can be written in 
the form (jl.2p and by setting e = it formally reduces to the Euler equation. For the time 
being, we focus on the Cauchy problem obtained by coupling (11.21) with the initial datum 

u(0, x) = uq(x). (1.3) 

Establishing general results on the convergence of the family of functions u e satisfying (11. 2J) , (11.31) 
is still a mayor open problem, but partial results have been achieved and ODE techniques 
have often played an important role: see again Dafermos [7] for a complete overview and a list 
of references. Let us just mention the work by Bianchini and Bressan [3], which established 
convergence (under suitable assumptions) of the family of functions u e solving (11. 2p . (11.31) in 
the case of the "artificial viscosity" B(u E ) = I (identity matrix). The proof in [3] employs a 
special decomposition of the gradient u £ x obtained by relying on ODE analysis. 

These notes aim at providing an informal overview on some of the ODE techniques that 
are more frequently applied to the study of the limit e — > + of (11.21) . The exposition is 
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addressed to non-experts, so notions are usually first introduced in a simplified context and 
then discussed in more general situations. Also, in most cases these notes only provide 
an heuristic idea of the proof of the the results that are introduced, and refer to books or 
to original research papers for a more rigorous discussion. The discussion has, of course, 
no sake of completeness: in particular, to simplify the exposition very few references are 
provided. For a more satisfactory bibliography, one can refer to the books by Dafermos [7J 
and by Serre [19] (conservation laws and their viscous approximation). An extremely rich 
exposition of the analysis of invariant manifolds for ODEs is in the book by Katok and 
Hasselblatt [H] (see also the book by Perko [17]). 

To simplify the exposition, most of the following sections (all but the last one) focus on 
the case of the "artificial viscosity" B(u e ) = /, so that (11.21) reduces to 



However, many of the considerations extend to more general cases (sometimes the extension 
is not straightforward, though). The details concerning this extension can be found in the 
original research papers that are quoted in the following sections. 

Also, most of the analysis described in the following sections apply to the study of the 
non-conservative case 



where A is a suitable N x N matrices. Note that (11.2j) can be written in the form (II. 5p by 
setting A(u £ ,eu s x ) = Df(u e ) — eB(u £ ) x , where Df denotes the Jacobian matrix of the flux 
/. Again, the details concerning this extension are available in the original research papers 
quoted in the following sections. 

The exposition is organized as follows: Section [2] informally describes the connection 
between the analysis of the limit e — > + of the viscous approximation (11.21) and the study 
of ODEs by introducing the notions of traveling waves and boundary layers. Section [3] 
deals with the Center Manifold Theorem and Section [4] mentions some applications to the 
study of the viscous approximation. Section deals with the Stable Manifold Theorem, 
while the last section of these notes (Section E]) is more technical and describes a possible 
way of extending the definition of center and stable manifold to a class of singular ordinary 
differential equations arising in the study of the viscous approximation defined by the Navier- 
Stokes equation in one space variable. 

Acknowledgments. These notes were originally prepared for a course given at the Seventh 
Meeting on Hyperbolic Conservation Laws and Fluid Dynamics: Recent Results and Research 
Perspectives. The conference was held at SISSA, Trieste and was organized by Fabio Ancona, 
Stefano Bianchini, Rinaldo M. Colombo and Andrea Marson. The author warmly thank them 
for the kind invitation. 



u\ + f{u e ) x = eu% 



(1.4) 



u\ + A(u £ , eu £ x )u £ x = eB{u £ )u £ x 



(1.5) 
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2. Some motivations 

This section aims at discussing some links between the study of the viscous approximation 
of a system of conservation laws in one space dimension and the ODE analysis by informally 
illustrating the notions of traveling waves and boundary layers. 

2.1. Traveling waves. Let us start with some heuristic considerations. For every e > 0, 
f ll.4j) is a parabolic equation and one expects a strong regularizing effect. Conversely, it 
known that, even if the initial datum uq is smooth, nevertheless in general a classical solution 
(i.e., a solution of class C 1 ) of the Cauchy problem (11.11) . (11.31) may break down in finite time: 
remarkably, this may happen even in the scalar case N = 1 (see e.g. Dafermos [TJ Chapter 
4.2]). Hence, when studying the limit e — > + one handles a family of very regular functions 
which is expected to converge to something which has a much weaker regularity. Traveling 
waves are often regarded as useful tools to handle this behavior and to describe the formation 
of singularities in the limit. 

To fix the notations, let us introduce the following definition. 

Definition 2.1. Let u + ,vT E R N and a E R be given. The function U : R -> R N is a 
traveling wave for (11.21) joining the states u + ,u~ E R n and having speed a if U(y) satisfies 

U" =[f{U)-aU]' (2.1) 

and 

lim = u~ lim = u + (2.2) 

y— *— oo y~ »+oo 

In the previous expression, U' and U" denote respectively the first and the second deriv- 
ative of U with respect to the variable y. Some observations are here in order. First, the 
link between (12. ip and (ll.4p is the following: if U solves (12. ip . then one can verify by direct 
check that a solution of (11.4j) is obtained by setting 

u £ (t,x) ;=u(^—^j . (2.3) 

Also, the function U is a solution of the ordinary differential equation (12.11) . hence it is of 
class C 2 and the same regularity is inherited by u e . However, by taking the pointwise limit 
e — > + of u e we get 



u(t,x) := lim u £ (t,x) = lim U 

E->0+ 6^0+ 



x — at 



lim^^oo U(y) if x — at < 



limj / _ >+00 U(y) if x — at > 

One can then verify that the pointwise limit u(t, x) provides a distributional solution of (II. ip . 

Summing up, we have that, if conditions (12.21) hold, then the solution of the parabolic 
equation defined by (12.31) is regular for every e > 0, but converges pointwise to a discontinuous 
distributional solution of the conservation law. This is why traveling waves are often used 
to study the non trivial behavior mentioned above, namely the appearance of singularities 
when passing from the vanishing viscosity approximation to the hyperbolic limit. 
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Let us now why the ODE analysis comes into play: the following considerations are quite 
informal, for a more precise and rigorous exposition see for example Bianchini and Bressan [21 
Sections 3, 4]. Let us consider the equation ( 12. ip : usually, the values u + , u~ and a are not 
prescribed. Conversely, one wants to understand for which values of u + , u~ and a there is 
a solution of (12. ip satisfying 

lim = u~ lim = u + . 

y— »— oo y—i+oo 

To select a unique solution of (12. ip . one needs to assign 2iV + 1 parameters: for example, one 
can assign a Cauchy datum for U and for U' and the value of a (remember that U takes values 
in R ). As a matter of fact, in many situations one does not manage to assign so many 
conditions and apparently one ends up with an underdetermined problem. By assigning 
2N + 1 parameters, however, one is neglecting an important information: the values of u + 
and u~ are not prescribed, but nevertheless both the limits lim x ^± U (y) exist and are finite. 
In particular, this implies that U is bounded on the whole M. In many situations, one 
actually has even more information: given a value u G M. N , one looks for values of w + , u~ 
and a satisfying the following property. 

Both u + and u~ are close to u, a is close to Aj(w). (2.4) 

Here Aj(tt) denotes the i-th eigenvalue of the Jacobian matrix Df(u). As we see in Section[3l 
the Center Manifold Theorem allows to exploit the information (12.41) to reduce the number 
of conditions one needs to impose on (12.11) to select a unique solution. 

2.2. Boundary layers. Boundary layer phenomena are observed when studying the viscous 
approximation of initial boundary value problems. 

To illustrate the heart of the matter, let us focus on a specific situation. Consider the 
family of problems 

f u% + f{u e ) x = eu £ xx 

\ u £ (t,0) = u b (t) u £ (0,x) = u (x). 

As established in Ancona and Bianchini [1] , under suitable assumptions on the data uj, and 
Mo, for every t > the family u e (t,-) converges in L\ oc to u(t, •), a distributional solution 
of (II. ip satisfying the following condition: for every t > 0, ToWar^-u^, •) < +oo. The point 
where boundary layers come into play is the fact that, in general, the boundary condition is 
lost, namely (note that the limit below exists since the total variation is bounded) 

lim u(t,x) ^ u b (t). (2.5) 

x— »0+ 

See also the previous work by Gisclon [TD] for the proof of the local in time convergence for 
more general viscous approximations (11.21) . 

To study the "loss" of boundary condition (I2.5p . the key point is the analysis of the steady 
solutions of (11.41) . satisfying f{u) x = eu xx . To understand why steady solutions are able to 
capture this behavior (the loss of boundary condition), let us start by considering a toy 
model. 



VISCOUS APPROXIMATION VIA ODE ANALYSIS 



5 



Example 2.2. Consider the vanishing viscosity approximation of a scalar, linear conservation 
law 

u\ + au e x = eu e xx a < (2.6) 
and let us focus on the steady solutions u £ (t, x) = U £ (x), obtained by solving the linear ODE 

all' = eU'J. (2.7) 
Also, let us impose on (12.71) the conditions 

U £ (0) = u b lim U £ (x) = u . 

x— >+oo 

Then by performing explicit computations we get U £ (x) = [u b — Uo\e ax ^ £ + u . Hence, by 
letting e — > + and, by exploiting a < 0, we obtain that, for every x > 0, U £ (x) converges to 
U(x) = u and hence the boundary condition U £ (x) = u b is "lost" since linx,,^ + U(x) = u . 

Some observations are here in order: first, by introducing the change of variables (t, x) h- > 
(et,ex) one can transform the equation f(u £ ) x = eu e xx into f(u) x = u xx . Also, assume that 
we want to study the "loss of boundary condition" for a linear system of conservation laws 

Ut + Au x = u xx where u G M. N . 

Here, A is an iV x N matrix and we assume that it admits N real and distinct eigenvalues. 
We want to proceed as in Example 12.21 by focusing on the steady solutions u(x,t) = U(x) 
and assigning the conditions U(0) = Uf, and lim x ^ +00 U(x) = u . However, it turns out that 
these conditions are compatible if and only if the vector u b — u belongs to the subspace of 
WL N generated by the eigenvectors of A associated to strictly negative eigenvalues. Let us 
now consider the general case: the second order approximation of a non linear system of 
conservation laws 

Ut + f(u) x = u xx where u e W. N . 
To extend the argument in Example 12.21 , we have to study the system 

[f(U)]' = U" 

(2.8) 

U (0) = u b lim^+oo U (x) = w - 

In view of the previous considerations it is thus natural to ask ourself the question: what 
are the values of Uf, and uq ensuring that (12.81) admits a solution? Note that such a query 
has some similarities with the problem we discussed at the end of Section 12.11 namely we 
are concerned with an ODE, we want to prescribe some asymptotic behavior and we ask 
ourselves what are the data compatible with such an asymptotic behavior. 

As we will see, an answer to the above question can be obtained by relying on the Stable 
Manifold Theorem, discussed in Section 

Finally, to fix again the notations, we introduce the following definition. 

Definition 2.3. Let the states Ub, u e M. N be given. A boundary layer for (2.1) connecting 
the states it& and u is a solution of system (I2.8p . 
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3. The Center manifold Theorem 

In this section we are concerned with studying the asymptotic behavior of the solutions 
of the ordinary differential equation 

V' = G(V) where V eR d 

in a neighborhood of an equilibrium point V. Without any loss of generality we can assume 
V = 0, namely G(0) = 0. In Section H] we will then apply our analysis to the study of the 
system (12.11) . 

Let us start by considering two simple examples where G is linear. 
Example 3.1. Assume that V = (x,y,w,z) and that 
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(3.1) 



By solving the ODE explicitly we obtain 

/ x(t) \ j e 2t x(0) \ 

y(t) = e-*j/(0) 

w (t) w(0) cos 3t — z(0) sin 3t 

\ z(t) J \ z(0)cos3t + w;(0)sin3t / 

Note that the eigenvalues of the above matrix are 2,-1 and ±3z. If the initial datum belongs 
to the eigenspace associated to the eigenvalue with positive real part, namely if the initial 
datum is in the form 

(x(0),0,0,0), 

then the solution blows up exponentially fast when x — > +oo, while when x — > — oo it 
converges with exponential speed to the equilibrium point. Conversely, assume that the 
initial datum belongs to the eigenspace associated to the eigenvalue with negative real part, 
namely the initial datum is in the form 

(0,y(0),0,0). 

Then the solution converges exponentially fast to the equilibrium for x — > +oo, while it 
blows up exponentially fast when x — > — oo. Finally, assume that the initial datum belongs 
to the eigenspace associated with the eigenvalues with zero real part, namely it is in the 
form 

(0, 0,«;(0), *(())). 

Then the solution does not converge exponentially neither at +oo nor at — oo. However, the 
solution is globally bounded on the whole real line. 

The second example focuses on a linear system where, however, the matrix is not diago- 
nalizable. 
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Example 3.2. Assume that V = (x,y) and that 

: ; H-)a) 

The above matrix has eigenvalue with multiplicity 2 and hence the whole R 2 is the 
eigenspace associated with eigenvalues with zero real part. By solving the equation explicitly 
we get 

x \ _ f z(0) +y(0)t 

y J ~ V y(o) 

In this case the solution is not bounded on the whole R unless y(0) = 0. However, it does 
not blow up exponentially neither at x — > +oo nor at x — > —oo. 

3.1. The center manifold theorem. 

3.1.1. Statement of the theorem. Let us first introduce some notations: B(0,5) is the ball 
with radius 5 and center at in M. d . Also, DG(0) is the Jacobian matrix of G computed at 
V = and V c will be its center subspace, namely 

V c = Z 1 ®Z 2 ® ...Z nc , (3.3) 

where each of the Zi is an eigenspace associated to an eigenvalue Aj with zero real part, 

Zi = {v such that [DG(6) - X^fv = for some k < d) ReA, = 0. 

We are now ready to introduce the main result of this section: 

Theorem 3.3. Consider the first order ODE 

V' = G(V), V G R d , (3.4) 

where G : IR d — > M d is a C 2 function satisfying G(0) = 0. Assume that the center space 
V c defined by ( 13.31) is non trivial, namely V c ^ {0}. Then there exists a constant 5 > 
small enough and a continuously differentiable center manifold A4 C satisfying the following 
conditions: 

(1) A4 C is parameterized by a function 

C : B(6,5) f]V c ^R d . 

Also, A4 C is tangent to V c at the origin. 

(2) Ai c is locally invariant for the ODE (13.41) . Namely, ifVo £ M. c , then the solution of 
the Cauchy problem 

' V = G(V) 
V(0) = V 

belongs to Ai c if \t\ is small enough. 

(3) let V(t) be a solution of fl33J) such that V(t) E B(5,6) for every t E R. Then 

V(t) E M c for every t E R. 
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An important remark is that the center manifold (I3.4p about the equilibrium point is not 
unique (see for example the lecture notes by Bressan [6J for an explicit example). In other 
words, suitable examples show that in general there is more than one manifold satisfying 
conditions (1)..(4) in the statement of the theorem. We will come back to this point in the 
proof of the theorem, when it will be clear where this lack of uniqueness come from. 

The proof of Theorem 13.31 discussed in here is the same as in the notes by Bressan [5J. For 
a more general viewpoint, one can refer to the book by Katok and Hasselblatt [T4"] . 

3.2. Proof of Theorem 13.31 The proof is divided in several steps, and only some of them 
are provided below (see [6J for the complete argument). 

Step 1: heuristic. The heuristic idea is that the manifold M c should play the role that 
in the linear case is played by the center space. In exploiting this idea one should keep in 
mind that, as Example 13.21 shows, the solutions of a linear system lying on the center space 
in general are not bounded on the whole real line: the only requirement we can reasonably 
impose is that they did not blow up exponentially fast neither at +oo nor at — oo. 

To make an extension to the general non linear case we have to handle some technical 
difficulties. First of all, we need to introduce a localization argument. Also, because of the 
non linearity we cannot require that the solutions lying on the center manifold do not blow 
exponentially neither at +oo nor at — oo, but only that they are controlled by exp(rj\t\), 
where rj is a small positive constant. 
Introduction of a cut-off function Let 

p(V) = G{V) - DG(0)V 

be the second order term in the expansion of G. We then have 

\p(V)\<C\V\ 2 and \Dp(V)\ < CV (3.5) 

for a suitable constant C and for V in a small enough neighborhood of 0. 

We want to reduce to the case p has compact support and small C 1 norm. We thus 
introduce a smooth, even cut-off function p : K. — > R satisfying 

< p(t) < 1 for every t e K and p(t) = 

We then define 

ps(V) :=p(V)p(\V\/5), 

thus obtaining that the support of p$ is contained in B(0,28). In the previous expression, 
5 denotes a small constant whose exact value is determined in the following. By relying 
on (13.51) we also get 

IMIc°< sup \p(V)\<AC5 2 

XeB{6,28) 

and 

\Dp s (V)\ = \Dp(V)\p+\p(V)\\p'\ 1 -<C5. 

o 





if \t\ 


o 


if \t\ 



< 1 

> 2 
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By combing the previous observations we can make the C 1 norm arbitrarily small by choosing 
5 small enough. 

In the following, we will study the solutions of the ODE 

V = DG(0)V + Ps (V). (3.6) 

By construction, 

DG(0)V + p s (V) = G(V) when \V\ < 5 

and hence any claim holding for the solutions of (13.61) holds for the solutions of the original 
equation V = G(V) confined in B(0,5). 

Step 3: definition of a fixed point problem We first have to introduce some further 
notations: for simplicity, in the following we denote by A the Jacobian matrix DG(0). Also, 
the stable and unstable space of A are the subspaces of WL d associated to the eigenvalues with 
negative and positive real part respectively. Namely, 

V s = Nj. © N 2 © . . . N ns V u = P 1 © P 2 © . . . P nu (3.7) 

where the iVj and the Pj are eigenspaces associated with eigenvalues with negative and 
positive real part respectively: 

N { = {v G R d : (A - \l) k v = for some k < d} ReAj < 

and 

Pi = {v E R d : (A - Xilfv = for some k < d} ReXi > 0. 
In particular, since 

R d = V s © V u © V c , 

then any x G lR d satisfies x = tt s x + tt u x + ir c x, where 

tt s : R d -> V s n u : R d -> V u ir c : R d -> V c 

are the projections onto V s , V u and respectively. 
Note that the ODE 

V = AV+p s (V) 
can be then written as the coupling of the ODEs 

' (ir a V)' = A7r s V + 7T sP5 (V) 
< (n u V)' = An u V + 7T uP5 {V) 

k (tTc^)' = AlX c V + 7T c p5(V). 

We have exploited the following equalities: 

tt s AV = Att s V tt u AV = Att u V tt c AV = Att c V for every V G R d 
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We solve each equation separately and by applying the "variation of constants" formula we 
obtain that the components ir s V(t), TT u V{t) and TT c V(t) satisfy 

7T s V(t) = e A ^7r s V(t s ) + f e A ^7T s p 5 (V(s))ds 



n u V(t) 



7l c V(t) 



,A(t-s) 



Tr u ps(V(s))ds 



J t u 

z A ^v c V(t c ) + f e A ^7r cPs {V(s))ds. 



(3.8) 



In the previous expression, t s , t u and t c denote real values which will be assigned in the 
following. Also, we have exploited the standard notation 



exp{At) = 



t n A r 



n=0 



Let Aj denote as before an eigenvalue of A, then we set 

/?+ := min{ReAj : ReAj > 0} 

and 

(3_ := min{-ReAj : ReAj < 0}. 
There exists some constant C > such that, for every v G lR d , we have 

\e At 7T u v\ < Ce p+t,2 \v\ for every t < and \e At 7r s v\ < Ce~ p - t/2 \v\ for every t > 0. (3.9) 

We now fix a constant positive constant rj > satisfying 

■p + (3. 



rj < mm 



and set 



Y v := {V G C°(] - oo, +oo[; R d ] such that sup \V{t)\e~^ < +oo}. 
Then Y v equipped with the norm 

||V||„ := supple"" 1 ' 1 
tent 

is a Banach space. Note that 

\V(t)\ < e^H^H, for every t G R, V G Y„. 

Loosely speaking, the space Y n contains all the function which "do not blow up with too fast 
exponential speed" neither at +oo nor at — oo. 

The goal is now defining the map C by solving a suitable fixed point problem in Y v . First, 
we observe that, if the function V G Y v , then by relying on (13.91) we get that, for any t s < t, 

\e A{t -^7i s V(t s )\ < e-^^e-^'WVWr, = e-^ 2 \\V\\ v e^^-^ 
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and hence 

lim \e A{t ~ ts) ir s V(t s ) \ = for any given tel. 
In the same way one can show that, if the function V G Y v , then 

lim \e A{t ~ tu) n u V{t u )\ = for any given tel. 

We now go back to formula (13.81) . we set t c = 0, x c = 7r c V(0) and we let t s — > — oo and 
t u — > +oo, eventually obtaining 

V(t) =TT c V(t) + TT s V(t) + 7T u V(t) = e At x c + f e Ait ~^7r cP5 (V(s))ds 







t rt 



(3.10) 



+ / e A ^ s \ s p s {V{s))ds + / e A ^n uPs (V(s))ds 

J — oo J +oo 

for every V e Y v . 

Further steps in the proof of the Center Manifold Theorem 

• By applying the Contraction Map Theorem, one can show that, for any x c e V c , 
the fixed point problem (13.101) admits a unique solution V G Y v . To verify that the 
hypotheses of the Contraction Map Theorem are satisfied one has to exploit that p$ 
is a function with compact support and small C l norm. 

• The map 

C : V c -> M d 

which parameterizes the center manifold is then defined by setting 

<f> c {x c ):=V{0). (3.11) 

• One can then prove that such a map is continuously differentiable and that the 
tangent space at the origin is V c . 

• By construction, the map (j) c satisfies the following properties: first, for any x c G V c , 
the solution of the Cauchy problem 

V' = AV+p s (V) 

(3.12) 

V(0) = 4>c{x c ) 

belongs to Y v , and hence it "does not blow up with fast exponential speed" neither 
at +oo nor at — oo. Conversely, any solution of A3. 12[) verifying 

supll/^le"" 1 * 1 < +oo 
t 

satisfies V(0) = (j> c {x c ) for some x c G V c . 

We set A4 C := 4> C {V C ) and we eventually come back to the original equation V = 
G(V). As pointed out before, G(V) = DG(0)V + p 5 (V) when \ V\ < 5. Consequently, 
any claim concerning the solution of V = DG(0)V + ps(V) can be reformulated in a 
claim concerning the solutions of V = G(V) which are confined in -8(0, 5) for every 
tel. By relying on this observation one concludes the proof of Theorem 13. 3[ 



12 



LAURA V. SPINOLO 



A remark concerning non-uniqueness As pointed out before, in general the center 
manifold of (13.41) about a given equilibrium point is not unique. The reason is that in the 
proof we have some freedom in choosing the cut-off function p: hence, both the map 4> c 
defined as in (I3.1ip and the manifold A4 C depend on the choice of p and are not unique. 

4. Some applications of the Center Manifold Theorem to the analysis of 

traveling waves 

The Center Manifold Theorem has been applied by several authors to the analysis travel- 
ing waves: see Dafermos [7] for a list of references. In this section we will discuss a possible 
approach, due to Bianchini and Bressan [3] and in the second part we will deal with appli- 
cations to the study of the vanishing viscosity approximation u\ + f(u £ ) x = su £ xx . 

4.1. A center manifold of traveling waves. Let us now go back to the problem we 
discussed at the end of Section 12.11 we are given the traveling wave equation 

U" = [f(U) - all]' where U G R N (4.1) 

and a value u G M. N . We want find solutions U admitting finite limits both at x — > ±oo, 
lim^-i-oo = m ± and satisfying the following properties: first, for every y G R, U(y) belongs 
to a small neighborhood of u. Also, a is close to \{u), where \{u) is a given eigenvalue of 
the Jacobian matrix Df(u). More precisely, we want to exploit the information we have on 
U to reduce the number of conditions one has to assign on (14.11) to select a unique solution. 
By setting p := U' we get that (14. ip can be written as the first order system 

U' = p 

p'= [Df{U)-al]p (4.2) 
a' = 0. 

We now apply the Center Manifold Theorem: first, we observe that the point V = [U, 0, \i(U)) 
is an equilibrium for (14.21) . Hence, from the third claim in the statement of Theorem 13.31 we 
deduce the following. There exists a constant 5 > and a center manifold Ai c such that 
any solution of (14.21) satisfying 

| (U(y) , p(y) , a) - V | < 5 for every y G R (4.3) 

lies on Ai c . As we mentioned in Section [21 the center manifold of a system about a given 
equilibrium point is in general non-unique. In the following we will arbitrarily fix one center 
manifold and denote it by A4 C . 

We now want to investigate the structure of Ai c . First, we compute its dimension, which 
is the dimension of the center space. By setting 



V=(U,p,aY and G{V) = (p, [Df(U)~ al]p, 0)* 
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we get that V G M 27V+1 and the Jacobian matrix DG(V) is given by 

/0 I 0\ 

DG(V) = Df(u) - \{u)I • (4.4) 
V / 

In the previous expression, and / denote the null and the identity NxN matrix respectively. 
Let us assume that the matrix Df(u) has N real and distinct eigenvalues: this hypothesis 
of so-called strict hyperbolicity is quite standard in the study of conservation laws (see for 
example Dafermos Section 9.6] for a discussion on the "bad behaviors" exhibited by 
systems violating strict hyperbolicity). Let f denote a unit eigenvector associated to Xi(u), 
namely 

\Df(u) — \i{u)\fi = and \fi\ = 1. 
Then the center space of the matrix (14.41) is 

V c = j(u,p, a) such that u G M. N , p = af for some aGl, and cr G M j. 

Hence, the dimension of V c and of Ai c is N + 2. This gives an answer to one of the questions 
we asked ourselves: if we are looking for solutions of (14.21) satisfying (14.31) . we need to impose 
N + 2 conditions (instead of 2N + 1) to select a unique solution. 

Let us get some more precise information about the structure of Ai c . Let 

<f) c : B{6,5)nV c ^R 2N+1 (4.5) 

be a parameterization of Ai c . From the proof of Theorem 13.31 we infer that <p c can be chosen 
in such a way that n c o(j) c (x c ) = x c for any x c G V c . By applying this property to system (14.21) 
we get 

tt c o (j) c (U, afi, a) = (U, af h a) for every U G R N , a G R, a G R. (4.6) 

and hence 

4> C (U, af, cr) = (V, ?/> c ([/, a, cr),o), 

where ip c : V c — > R N is a suitable function whose properties we want to investigate. Let 
{r*! . . . r/v} be a basis of R N entirely made by unit eigenvectors of Df(u). Then we denote 
by ipi c . . .ipNc the corresponding components of %fj c and we get 

N 

i[) c (U, a, cr) = ^ ^kc{U, a, o)f k . (4.7) 

k=l 

To write (14. 7p in a more convenient form, we first exploit (14. 6 p to get ipi C {U, a, a) = a. Also, 
any point (U, 0, a) is an equilibrium for (I4.2p . so it belongs to A4 C . This implies that 

TV 

= i'kciU, a, a)fk and hence that ipkc(U, 0, cr) = for every k, U and a. 
k=i 
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By exploiting the regularity of the map ipkc, we deduce that 4>kc(U, a, a) = a^jU, a, a) for 
a suitable function ipkc- By combining the previous observations we then get that (14.71) can 
be rewritten as 

p = a (?i + ^2 i>kc(U, a, a)f fc ) . 

k^i 

Since Ai c is tangent at the equilibrium point to the center space V c we also have 

( y r i + J2^kc(U,0,X l (U))f k }j =f h (4.8) 

namely ipkc(U, 0, \i(U)) = for every k ^ i. By setting 



"i(U, a, a) := (f, + ^ 4>kc{U, a, a) 



Ti{U, a, a) ■= \ n + Vkc\U, a, °Fk 

k^i 

we then obtain that (14. 7[) can be rewritten as 

p = afi(U,a,a). 

For convenience, we want to reduce to the case of unit vectors, so we write the above relation 
as 

p = a\ri(U,a,a)\ -— rr, 

\ri(U,a,a)\ 

where \fi(U, a, a) | denotes the length of the vector. Note that (14.81) implies that fi(U, 0, Xi{U)) 
fi. Hence, the map 

a i — ^ v(a) = a\fi(U, a, a)\ 

is locally invertible in a neighborhood of (U, 0, Xi(U)): by assuming that the size of the 
constant 5 in (I4.5[) is small enough, we eventually get that the manifold M. c is described by 
the relation 

p = vf-i(U,v,a), (4.9) 

where we have set fj = fj/|fj|. 

4.2. Center manifold and convergence of the vanishing viscosity approximation. 

In this section we just mention the fact that center manifold techniques have been exploited 
by Bianchini and Bressan [3] to show the convergence of the solutions of the family of Cauchy 
problems 



xx 



.«■«>,*) = «o(*). (4 - 10) 
More precisely, the following holds. 

Theorem 4.1 (Theorem 1 page 229 in [5]). Assume that the flux function f : WL N — > M. N is 

of class C 2 and that, for every u G M. N , the jacobian matrix Df(u) has N real and distinct 
eigenvalues. There exists a small enough constant 5 > such that, if the initial datum 
in (14.101) satisfies 

u G L x (R) and TotVar{w } < 5, 
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then the following properties hold. 

(1) For any e > 0, the Cauchy problem (14.101) admits a solution defined for any t > 0. 
Also, there exists a constant C, which does not depend neither on e nor on t, such 
that 

TotVar^w^t, •)} < C5. (4.11) 

(2) The solution u s is stable in L 1 with respect to time and to L 1 perturbation of the 
initial data. Namely, if by using a semigroup notation we denote by t t— > Sfuo the 
solution of (14.101) . then 



\\S^Uq — S^oIIl 1 < -^^ll^o — v \\ L i + \t — s\ + \ Vei — y/es]) for every t, s > 0. 

In the previous expression, L denotes a constant which does not depend neither on e 
nor on the time. 

(3) As e — > + , for every t > 0, S^uo converges strongly in Lj oc (R) to a function S t uo 
providing a distributional solution of (14.101) . Also, the semigroup t \— > S t uo enjoys 
the stability property 

\\S t u — S s v \\ L i < L^\\u — v \\ L i + \t — s\j for every t, s > 0. 

4.3. Solutions of the Riemann problem. In this section, as a direct application of the 
center manifold construction discussed in Section 14.1} we go over the analysis of the solution 
of the so-called Riemann problem 



u t + f(u) x = 

u + x > (4.12) 
u~ x < 



u(0, x) 



Here u + and u— are two given values in M, N and in the following we focus on the case 
\u + — u~\ < 5 for a small enough constant 5. The Riemann problem is a very specific example 
of Cauchy problem, but is nonetheless extremely important. Indeed, its study provides the 
building block for the construction of approximation schemes (Glimm scheme |TTJ , wave 
front-tracking algorithm) that have been used to establish global existence and uniqueness 
results for the system of conservation laws u t + f(u) x = 0: see for example Dafermos [TJ 
Chapter 14] and the references therein. 

The first issue one has to address when studying (14.121) is that, in general, a distribu- 
tional solution of (14.121) is not unique: in the attempt to select a unique solution, various 
admissibility criteria have been introduced. An admissibility criterium was introduced by 
Lax [IS], who also constructed an admissible solution of (14.121) . The analysis in (TS] relied 
on some technical assumptions on the flux / which were later relaxed in Liu [16]. See also 
Tzavaras [20] and Joseph and LeFloch [13]. Bianchini and Bressan [3], Section 14] provided 
a general construction of an admissible solution obtained by taking the limit e — > + of the 
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vanishing viscosity approximation 




eU xx 



u + x > 
u~ x < 

In both [151 dSj and [3] the key point is the construction of the so-called i-wave fan curve. 
Loosely speaking, the basic idea is the following: given u~ G M. N , the i-wave fan curve 
Ti(u~, s) contains all the states u + such that an admissible solution of the Riemann prob- 
lem ( 14.121) is obtained by gluing shocks (or contact discontinuities) and rarefaction waves 
belonging to the z-family. In particular, this implies that there is a self-similar, admissible 
solution u(t, x) = U(x/t) of (I4.12p such that U(£) varies only in a small neighborhood of the 
eigenvalue \i(u~) and is constant elsewhere, U(£) = u~ for £ < \i(u~) — 5 and U(£) = u + 
for £ > \i{u ) + 5. 

We now informally go over the main steps of the construction in [3]. 

(1) We fix u~ G ~R N and we consider the same fj as in (14. 9p . obtained by applying the 
Center Manifold Theorem to system (j4.2p about the equilibrium point (u~ , 0, Aj(u~)) . 
Here Aj(u~) is the z-eigenvalue of Df(u~). Given s > 0, we consider the following 
fixed point problem, defined on the interval [0, s]: 

u(t)=U~+ fi{u(g),Vi(g),<Ti(g))d£ 



Vi(r) = /i(r, it, Ui, <7i) - conv/i(r) (4.13) 

o-i(r) = — conv/j(r). 
ctr 

Here 

fi( T )= K( u (C),Vi(C),^i(C))dC and A*(u(£), ^(0) = n), 

namely Aj can be regarded as a generalized eigenvalue. Also, conv/j(r) denotes the 
convex envelope of fi(r), 

conv/j(r) = sup |(?(t) s.t. g is convex, < /(£) for every £ G [0, s] j. 

By relying on the Contraction Map Theorem, one can show that (14.13p admits a 
unique continuous solution (u(t), Uj(r), ct^t)) defined in a small enough neighbor- 
hood of (u~, 0, Aj(w~). The functions (u(r), Vi(r), <Ti(r)) are all defined on the interval 

(2) Assume that the interval [a, b] C [0, s] satisfies the following property: for every 
r G [a, 6], fi(r) = conv/j(r). Then one can show that fj(u(r), Ui(r), 0i(t)) = rj(u(r)) 
for every r G [a, 6] and moreover that the solution of the Riemann problem 

Ut + f(u) x = 

, , j u(b) x > 
u(0, x = < ) < „ 
v ' ' \ u(a) x < 
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is a rarefaction wave of the i-th family (hence, in particular, the solution is continu- 
ously different iable in ]0, +oo[xR). 
(3) Conversely, assume that for every r g]c, g?[C [0, s], fiij) > conv/^r) and that fi(c) = 
conv/i(c), fi(d) = conv/j((i). Then, in particular, one can show that Oi is constant 
on [c, d], Gi = Oi. Also, Vi{r) > on ]c, d[ and one can define the change of variables 

dx 1 

dt Vi(t) 

(n\ d ~ c 
x{0) - 



2 

mapping x :]c, d[— >] — oo, +oo[. Let r(x) denote its inverse: one can show that the 
function U(x) := u(t(x)) satisfies 

U" = (f(U) - OiU)' and lim U(x) = u(c), lim U(x) = u{d). 

X— > — oo x— > + oo 

In other words, there is a traveling wave connecting u(c) and and hence a solution 
of the Riemann problem 



Ut + f(u) x = 



u(0, x) 



u(d) x > 
w(c) x < 



is 

u(d) x > ait 
u(c) x < ad, 



u(t, x) 



which is called either shock or contact discontinuity depending on the values of 
Aj(u(<i)) and Aj(u(c)) (see for example the book by Dafermos [TJ Page 213] for the 
exact definition). 

(4) By relying on the previous steps and on further approximation arguments, one can 
show that if (u(r) , v^r) , crj(r)) is the solution of (14. 13j) . then the family u e solving 




£U XX 



u x < 
u(s) x > 



converges e — >• + to the function 



n~ if x < (Tj(0)t 
u(t,x)=^ w(r) if x = (jj(r)t 
w(s) if x > ai{s)t 

We then define the i-wave fan curve by setting Ti(u~, s) := u(s). 
(5) When s < 0, the construction of Ti(u~, s) is analogous, the only difference being that 
in ( 14.131) one has to take the concave envelope of fi instead of the convex one. 



18 



LAURA V. SPINOLO 



Finally, observe that the above construction can be extended to study the solution of the 
Riemann problem (14.121) obtained by taking the limit e — > + of the more general viscous 
approximation 

ut + f{u% = e(B{ U z)ul)^ 

under quite general assumptions on the matrix B (see Bianchini [2]). For extensions to the 
analysis of initial-boundary value problems see Bianchini and Spinolo [5]. 

5. Stable Manifold Theorems 

5.1. The Stable Manifold Theorem. Let us first recall some notations: given the ODE 

V = G(V), (5.1) 

let V be an equilibrium point, which without any loss of generality is assumed to be 0. As 
before, we denote by V s and V u the stable and the unstable space respectively, defined as 
in (13.71) . Also, as usual, DG(0) denotes the Jacobian matrix of G computed at 0. 

Theorem 5.1. Consider the first order ODE (15.11) where G : ~R d — > ~R d is a C 2 function 
satisfying G(0) = 0. If V s is non trivial, namely V s ^ {0}, then there exists a constant 5 > 
small enough such there exists a continuously differentiable stable manifold M. s satisfying the 
following conditions. 

(1) M. s is parameterized by a function 

S : B(0,5)nV s ->R d . 

Also, M. s is tangent to V s at the origin. 

(2) M. s is locally invariant for the ODE (13.41) . Namely, ifVo G Ai s , then the solution of 
the Cauchy problem 

' V = G(V) 



V(0) = V (5 ' 2) 
satisfies V(t) G Ai s for every \t\ small enough. 

(3) IfVo G M. s then the solution of the Cauchy problem ( 15. 21) satisfies 

lim \V(t)\e ct/2 = 0, (5.3) 

t— >+oo 

where 

c = min{— ReXi such that Aj is an eigenvalue of DG(0), Re\ < 0} (5-4) 

(4) Conversely, ifV(t) is a solution of (15. ip satisfying (15. 3p . then V{t) G M. s for every 
t large enough. 

The proof of Theorem 15.11 can be obtained by relying on an argument somehow similar 
to the one exploited in the proof of Theorem 13.31 a detailed proof can be found in the 
book by Perko [17] (see also the book by Katok and Hasselblatt [13] for a more general 
viewpoint). One should keep in mind, however, that the stable manifold (i.e., a manifold 
satisfying properties (1), (2), (3) and (4) in the statement of Theorem 15.11) is unique, while 
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the center manifold is not. As a matter of fact, one can also define a global stable manifold 

as 

M gs = (J {V(t) solving (JO) and satisfying V(0) G M 3 }. 

In this way, we obtain a manifold which is globally invariant for ( 15. II) . Namely, if Vq G Ai gs , 
then the solution V(t) of ( 15.21) satisfies V(t) G M 9S for every t. 

5.2. Applications to the analysis of the boundary layers. Let us now go back to the 
original problem described at the end of Section [2721 : given u , we want to describe the values 
of Ub such that problem 

( U' = f(U) - /(n ) 

1 U(0) = u b lim^+oo U (x) = u 

admits a solution. By combining parts (2) and (3) in the statement of the theorem we get 
that the previous problem admits a solution provided Ub belongs to the stable manifold A4 S , 
whose dimension is equal to the number of the eigenvalues of Dfiuo) with negative real part. 
If all the eigenvalues of the Jacobian matrix Df(uo) have nonzero real part, then as a matter 
of fact one can also prove the converse implication: if the above problem admits a solution, 
then u b G M s . 

5.3. The slaving manifold of a manifold of equilibria. Loosely speaking, the Stable 
Manifold Theorem we discussed in the previous section describes all the solutions of (15.11) 
that when t — > +oo converge exponentially fast to the equilibrium point 0. 

Let us consider a more general situation: assume that there exists a continuously differ- 
entiable manifold £ C M. d which contains and is entirely made by equilibria: 

V V G £, G(V) = 0. 

As a side remark, note that we are not assuming that £ is the set of equilibria, namely we 
do not rule out the possibility that there are equilibria that do not belong to £. In several 
situations, we are led to ask ourself the following question: is there any way to describe 
the solutions of (15. ip that converge to some point in £ ? Note that by applying the Stable 
Manifold Theorem we describe the solutions that converge to a given point, here we just 
require that the limit exists and that it lies somewhere on £. 

A positive answer to the previous question is provided by the theorem below. 

Theorem 5.2. Consider the first order ODE (15.11) where G : M> d —> M> d is a C 2 function 
satisfying G(0) = 0. Denote by E the space tangent to £ at the origin. Also, as before let 
V s be the stable space of G at the origin, defined as in (13 .7p : assume that V s is non trivial, 
namely V s ^ {0}. Finally, let the constant c be as in (15.41) . Then there exists a constant 
5 > small enough and a continuously differentiable manifold M. us satisfying the following 
conditions. 



20 



LAURA V. SPINOLO 



(1) M. us is parameterized by a function 



: B(0,S) n [V s ®E 



KV -LJ J — ► 

and it is tangent to V s © E at the origin. 

(2) M. us is locally invariant for the ODE f l3.4j) . Namely, if Uq G Ai us , then the solution 
of the Cauchy problem 

satisfies V(t) G M. us for every \t\ small enough. 

(3) If Vq G Ai us then there exists V^o G £ such that the solution of the Cauchy prob- 
lem (15. 5p satisfies 

lim - \4o|e ct/2 = 0. (5.6) 

t— >+oo 

(4) Conversely, ifV(t) is a solution of (15. ip satisfying (15.61) /or some |V^o| < 5, t/ien 
V(i) G /or every t /ar^e enough. 

The manifold can be regarded as a uniformly stable manifold because it contains the 
solutions of (15. ip that converge to some point in S with exponential speed uniformly bounded 
from below. As a matter of fact, the idea behind the definition of uniformly stable manifold 
is a particular case of the idea of slaving manifold. Loosely speaking, the slaving manifold 
relative to some manifold A contains all the solutions of (15.11) that when t — > +oo approach 
A with exponential speed: see again the book by Katok and Hasselblatt [H] . We came back 
to this in Section [6] by considering in Lemma 16.91 another class of slaving manifolds. 

Theorem 15.21 can be proven by relying on the Hadamard Perron Theorem: see Ancona 
and Bianchini [1] and the book by Katok and Hasselblatt |14j . 

6. Invariant manifolds for a class of singular ordinary differential 

equations 

This section focuses on a class of ordinary differential equations in the form 

d I=w/ {n (6 ' 1) 

where V G M d , F : lR d — >• M d is a smooth function and £ is a real valued, smooth function. 
The singularity of the equation comes from the fact that ( can attain the value 0. 

6.1. Motivations. Singular ODEs like (16.11) may arise in the analysis of the viscous approx- 
imation 

u t + f(u) x = (B(u)u x ) x 

when the N x N matrix B(u) does not have full rank, which is the case of most of the 
physically relevant examples. In particular, let us focus on the case of the compressible 
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Navier-Stokes equation in one space variable: 



Pt + {pv) a 







(pv) t + {^pv 2 + pj 



(6.2) 



pe + P YJt + 



-pv + pe + p 



k6 x + vv v x 



Here the unknowns p(t, x), v(t, x) and 9(t, x) are the density of the fluid, the velocity of the 
particles in the fluid and the absolute temperature respectively. The function p = p(p, 6) > 
is the pressure and satisfies p p > 0, while e is the internal energy. In the following we will 
focus on the case of a polytropic gas, so that e satisfies e = RO/^ — 1), where R is the 
universal gas constant and 7 > 1 is a constant specific of the gas. Finally, v(p) > and 
k(p) > denote the viscosity and the heat conduction coefficients respectively. 

Let us now compute the equations satisfied by the traveling waves and the steady solu- 
tions of (16. 2p : for simplicity, in the following we focus on steady solutions only, but the 
considerations we make can be repeated with small changes in the case of traveling waves. 
We have to consider the equation 



(pv) x = 
pv 2 +p) 

J x 



(6.3) 



-pv + pe + p 



k9 x + uvv 3 



To write (I6.3P in a convenient form we set 



w 



{v x , 9 X ) 



and hence we get 



v A 12 
A21 A22 



Px 
W 




b 



Pxx 
W x 



(6.4) 



Here, v4 12 G M 2 is a row vector, A 2 \ G M 2 is a column vector and they both depend on 
(p,v,9,w). Finally, A 22 and b are both 2x2 matrices depending on (p,v,6,w). The exact 
expression of these terms is not important here: the important thing to keep in mind is that 
the block b is invertible. 

To write ( 16.41) in a more explicit form, we first assume v 7^ and from the first line we get 



P.i 



— A 2 i ■ w, 



where • denotes the standard scalar product. By plugging this expression in the last two 
lines of (16 .4p we get 

b" 1 A 22 - - A 12 A 21 
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Note that this expression makes sense since b is invertible. By combining the previous 

considerations we get that (I6.3j) can be written in the form (16.11) provided that 

{ P\ 

v 



V 



-A 12 ■ w 
F(V) =| v w 

b^ 1 [A22V - A21A12] w 





\w ) 

and C{V) = v G R. Note that v is the velocity of the fluid and in general it can attain the 
value 0. 



Remark 6.1. Equation (16. 2p is the compressible Navier-Stokes equation written in Eulerian 
coordinates. If we write the Navier-Stokes equation using Lagrangian coordinates and we 
compute the equation satisfied by the traveling waves and by the steady solutions, we obtain 
standard ODEs with no singularity. See for example Rousset [18] for related analysis. 

6.2. What in principle can go wrong with (16. ip . There is a wide literature concerning 
the analysis of equations in the form 

% = -**{V), V G IR d , (6.5) 
at e 

where e is a parameter that can go to and F £ is a family of smooth functions. For an 
overview, see for example the lecture notes by Jones [12], which in particular provide an 
overview of Fenichel's works [H [9]. The reason why we cannot directly apply these results to 
the analysis of (16.11) is because the singularity ((V) depends on the solution V itself, while 
in ( 16.51) the singularity is just a parameter. Hence, when handling ( 16. II) we have to tackle 
the possibility that z(V) is nonzero at time t = 0, but then z(V) reaches the singular value 
in finite time. When this happens, the solution may experience a loss of regularity, as the 
example below shows. 

Example 6.2. Let us consider the ODE 

dv 1 /dt = -v 2 /v 1 , 66 , 
dv2/dt = —v 2 - ' 

It can be written in form (16. ip provided that V — (v±, v 2 Y, C(V) = v i an d 

F(V) = 



-v 2 

-V 2 Vi 



The solution of (16.61) is 



vi(t) = Jv 1 (0)+V2(0){e-t-l) 

V (6-7) 

v 2 (t) = v 2 (0)e~ t 

By choosing v 2 (0) > i'i(0) > 0, we get that C(^0 = v i(t) can attain the singular value for 
a finite t. Note that at that point t the first derivative dvi/dt blows up: thus, in particular, 
the solution (16.71) of (16. 6p is not C 1 . 
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6.3. Goals in studying ( 16. ip . The analysis concerning (16 .ip and discussed in the following 
aims at applying the ODE techniques described in the previous sections to the study of the 
viscous profiles of the Navier-Stokes equation. 

More precisely, the goals are the following. 

(1) Find possible extensions of the Center Manifold Theorem 13.31 and of the Uniformly 
Stable Manifold Theorem 15.21 to the case of the singular ODE (16. II) . Note that 
both theorems provide local results, so to study this extension we can restrict to the 
solutions of (16. ip belonging to a small enough neighborhood of an equilibrium point. 

(2) Rule out the losses of regularity discussed in Example 16. 2[ To do so, we have to find 
conditions that guarantee that the following property is satisfied. 

If £(y) ^ at t = 0, then ((V) ^ for every t > 0. (6.8) 

(3) Impose on ( 16.11) hypotheses satisfied by the viscous profiles of the Navier-Stokes 
equation. 

6.4. Hypotheses. This section lists the conditions imposed on the functions F and ( to 
study (16. ip . and comments on them. 

Hypothesis 1. The functions F : M. d — > M. d and £ : M d — > R are both regular (being of class 
C 2 is enough). Also, F{0) = and C(0) = 0. 

The second part of Hypothesis [1] simply says that we are assuming that our equation is 
singular and that the equilibrium point is 0, which of course is not restrictive. 

Hypothesis 2. We have V<C(0) ^ 0. 

Let S be the singular set 

S:={V: C(V)=0}. 

By relying on the Implicit Function Theorem and on Hypothesis [2] we get that in a small 
enough neighbourhood of the set S is actually a manifold of dimension d — 1, where d is 
the dimension of V. 

Hypothesis 3. Let Ai c be any center manifold for 

^ = F(V) (6.9) 
dr 

about the equilibrium point 0. If \V\ is sufficiently small and V belongs to the intersection 
A4 C PI S , then V is an equilibrium for (I6.9p . namely F(V) = . 

To understand why we impose Hypothesis [3] let us consider the following simple linear 
example. 

Example 6.3. 

dvi/dr = t> 2 /e 

dv 2 /dr = —Vi/s (6.10) 
dejdr = 0. 
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The first component of the solution is 

vi(t) = Acos(t/e) + Bsm(t/e), 

where A and B are real parameters. Letting e — > + , we get that in general if t ^ there 
exists no pointwise limit of v\. Note that (16.101) does not satisfy Hypothesis El Indeed, the 
center space is the whole R 3 and it coincides with the center manifold. However, 

F(V) = 

is not identically zero when e — 0. 

We now assume that there exists a non-degenerate set of equilibria. 
Hypothesis 4. There exists a one-dimensional manifold of equilibria 8 for the equation 

CLT 

which contains and which is transversal to S. 

The one-dimensional manifold £ is transversal to the manifold S if the intersection Sr\Ai eq 
has dimension zero. 

Hypothesis 5. For every V G S, 

V((V)-F(V) = 0. 

Hypothesis [5] is needed to rule out losses of regularity like the one in (16. 6p . Indeed, let us 
go back to example (16.61) : we have 

F(vx,v 2 ) = ( J , C(vi,v 2 ) = v 1 and S = {(ui, V2) : «i = 0}, 

hence 

VC(t'i,f2) • -F(wi,t; 2 ) = -v 2 , 
which in general is nonzero on S. 

To introduce the last hypothesis, let us consider the function 

q(y) .- vav)- F (v) 

- COO ' 

Thanks to Hypothesis [5] and to the regularity of the functions ( and F, the function G can 
be defined by continuity on the surface S. 

Hypothesis 6. Let U 6 5 be an equilibrium for (16.91) . namely C(V) = and F(V) = 0. Then 

G{V) = 0. 

Like Hypothesis [5], Hypothesis [6] is needed to rule out losses of derivatives. Indeed, let us 
consider the following example. 
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Example 6.4. Let us consider the ODE 



dvi/dt 
dv 2 I dt 
dv 3 /dt 



-v 3 

-v 2 /vi 
-v 3 , 



(6.11) 



then we can set 



v 2 , v 3 ) 




so that 



((vi, v 2 , v 3 ) = vi and S = {(vt, v 2 , v 3 ) : Vi = 0}, 



on S 



VC • F = -v 3Vl 

and hence Hypothesis [5]is satisfied. As a matter of fact, one can verify that Hypotheses [T]. . . H] 
are all satisfied as well. However, Hypothesis [6] is violated: indeed, 



G(v u v 2 , v 3 ) = -v 3 , 
which is in general different from on 

{(vi, v 2 , v 3 ) : such that F(vi, v 2 , v 3 ) = and v 2 , v 3 ) 
By solving (16.111) explicitly we get, after some computations, 

vi(t) =v 1 (0)(l-A)+Av 1 (0)e _t 



0} = {(0,0,^)}- 



v 2 {t) 



B 



1- A)e t + A 



i [(A-l)«i(0)] 



^i(0)e- 4 , 



where A and B are suitable constants depending on the initial data. If (A-l)vi(O) > 1, then 
the first derivative dv 2 /dt blows up at t = \n(A/A — 1). Note that this is exactly the value of 
t at which v±(t) attains 0. In general, for every fi(0) > if 1/(A — l)ui(0) is not a natural 
number, then the solution is not in C m for m = [1/{A — l)vi(0)} + 1. Here [1/(A — l)-ui(O)] 
denotes the entire part. Thus, we have a loss of regularity in higher derivatives. 

Remark 6.5. One can verify that Hypotheses [TJ . . [6] are all verified by the viscous profiles 
(travelling waves and steady solutions) of the Navier Stokes equation (16. 2p . 

6.5. A toy model. To get a flavor of the behaviors we expect to encounter, we start by 
considering the case when the singularity is only a parameter. Actually, to simplify things 
we focus on a toy model where everything is linear and we can carry out the computations 
explicitly. Let us consider the system 



dv i / dt 
dv 2 /dt 
de/dt = 



-5vi 
-v 2 /e 



(6.12) 



0. 
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First, we note that the subspace 



E = {(0, 0, e) : eel} 

is entirely made by equilibria. If e > the uniformly stable manifold relative to E is 
the whole space M 3 , since any solution of (16.121) decays exponentially fast to a point in E. 
However, 



-Bt 



and t> 2 



-t/e 



In other words, the speed of exponential decay of V\ is bounded with respect to e and v i is 
not affected by the presence of the singularity. Conversely, the speed of exponential decay of 
V2 gets faster and faster as e — * + . The second component of the solution is thus strongly 
affected by the presence of the singularity and can be regarded as a fast dynamic, while the 
first component is a slow dynamic. 

Summing up, in (I6.12p any orbit lying on the uniformly stable manifold relative to E 
decomposes as the sum of a slow and a fast dynamic. 

To see what happens if we look for extensions of the Center Manifold Theorem, let us 
modify (16. 12ft by adding two more equations: 



dvi/dt 
dv2 I dt 
dv^/dt 
dv^/dt 
de/dt = 



-5v i 

-v 2 /e 

-v 4 



(6.13) 



= ^3 

0. 



In this way, we have that 



v^{t),V4(t) ~ Acost + Bsint 

and hence we have a nontrivial center manifold. Note that we could not have added some- 
thing like 

dv^/dt = —Vi/e 
dv^/dt = v 3 /e 

because this would have violated Hypothesis [31 In other words, the center manifold can 
contain only slow dynamics. 



6.6. Slow and fast dynamics in the general case. Singling out the slow and the fast 
dynamics and defining the center and the uniformly stable manifold for the toy model (16.131) 
is a trivial task. However, to understand how we can proceed in the general case, let us first 
make some preliminary considerations, still focusing on (I6.13p . 

(1) We introduce the change of variables 



r = t/e. 
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Then (16 . 1 3[) becomes 

dv\/dr = —5viE 
dv 2 /dr = —v 2 

dv 3 /dr = —v^e (6.14) 
dv±/ dr = v 3 e 
de/dr = 0. 

(2) We can single out the slow dynamics of (16. 13ft by considering the center space 
of (16.141) . which is given by {(vi, 0, v 3 , v±, e)}. 

(3) Let us now restrict the analysis to the center space of (16.141) : we obtain 

dvi/dr = —5vie 
dv 3 /dr = —v^e 
dv^/dr = v 3 e 
de/dr = 0. 

Hence, we can come back to the original variable t, obtaining an equation with no 
singularity in it: 

dvi/dt = —5vi 

dv 3 /dt = -v± , Q15 , 
dvi/dt = v 3 \ ■ ) 

de/dt = 0. 

(4) We can then consider the center space of (I6.15p . which is given by V c = {(0, v 3 , v^, e)}. 
Then the center manifold of the original system (I6.13P is exactly V c . Also, we can 
observe that E = {(0, 0, 0, e)} is a subspace of equilibria for (I6.15P and we can consider 
the uniformly stable space of (I6.15P with respect to E: it is given by {(vi, 0, 0, s)}. 

(5) Let us now consider the uniformly stable manifold of (16. 13[) with respect to the 
manifold of equilibria E = {(0, 0, 0, e)}. It is given by {(i>i, v 2 , 0, 0, e)} : as mentioned 
before, any solution in there writes as a fast dynamic (in this case, (0,v 2 , 0,0)) plus 
a slow dynamic (in this case, (vi, 0, 0, e)) belonging to the uniformly stable manifold 
of (16.151) with respect to E. 

Let us now discuss how one can proceed in the general case 

dV 1 , s , 

* = W) F{n (616) 

Actually, we will not enter into the technical details, which can be found in pE]. The heuristic 
idea is proceeding by following the same steps (1). . . (5) as in the case of the toy model. Of 
course, now we have to handle additional difficulties coming from the fact that there is a 
non linearity in the function F(V). Also, ((V) is no more a parameter, but depends on the 
solution itself, and hence we have to rule out the possibility that ((V) ^ at t = but C(^0 
reaches in finite time. 
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Actually, as a preliminary step it is convenient to introduce a change of coordinates which 
allows to write (I6.16P in a more convenient form, however, this is a mostly technical step 
and in here we do not go into the details (they can be found in [IJ Proposition 4.1]). 

Let us now informally discuss how to extend the previous steps to the general case (16.16!) . 

(1) We would like to use the variable r defined by solving the Cauchy problem 



However, (16.171) defines a change of variables r :] — oo; +oo[^] — oo; +oo[ provided 
that ((V(t)) does not attain the value 0. Hence, we actually proceed as follows: we 
study the equation 



which is formally obtained from (16.161) by relying on (I6.17p . In particular, we focus 
on the solutions of (16.181) lying on suitable invariant manifolds. We then prove that 
for these solutions, the Cauchy problem 



defines a smooth diffeomorphism t :] — oo; +oo[^] — oo; +oo[. Hence, we can come 
back to the original equation (16.161) and a posteriori infer that the change of vari- 
ables (16.171) is well justified and it is inverse of the one defined by (16.191) . so our 
argument works. 

(2) Hypothesis [1] implies that F(0) = 0. We then consider the ODE (I6.18P linearized 
about the equilibrium point 0. We then say that a manifold of the slow dynamics is a 
center manifold of (I6.13P about the equilibrium point (any center manifold works). 

(3) By exploiting Hypothesis [61 one can show that if V(r) is a solution of (16.181) lying 
on the manifold of the slow dynamics, then the Cauchy problem (I6.19P defines a 
continuously differentiable diffeomorphism. Hence by restricting to the manifold of 
the slow dynamics we get that (I6.16P and (16.181) are equivalent. In particular, by 
relying on Hypothesis [3] we can show that the ODE satisfied by the solutions V(t) 
lying on the manifold of the slow dynamics is actually nonsingular. 

(4) We then consider system (I6.16P restricted on the manifold of the slow dynamics and 
we consider a center manifold about the equilibrium point 0: we can directly apply 
Theorem 13.31 because this system is nonsingular. This manifold can be regarded as 
a center manifold of the original system (16.161) : more precisely, we get the following 
theorem. 




(6.17) 



dV 
dr 



F(V) 



(6.18) 




(6.19) 



Theorem 6.6. Assume that Hypotheses^ ... are all satisfied. Then there exists a 
manifold A4 C defined in a small enough neighborhood ofO and satisfying the following 
conditions: 



VISCOUS APPROXIMATION VIA ODE ANALYSIS 



20 



(a) ifU(t) is a solution of ( I6TT6]) lying on M c and ((V) ^ at t = then ((V) ^ 
for every tel. 

(b) .M c s locally invariant for ( 16.1 6ft . Namely, if V G .M c ; i/ien £/ie solution of the 
Cauchy problem 



satisfies V(t) G .M c /or |t| small enough. 
(c) There exists a small enough constant 5 > s«c/j t/iai t/ie following holds. IfV(t) 
is a solution of (16.161) satisfying \ V(t)\ < 5 for every t G R ; t/ien \^(t) G .M c . 

Remark 6.7. In the statement of the Center Manifold Theorem 13.31 and of the Uni- 
formly Stable Manifold Theorem 15.21 one assumes that suitable subspaces (the center 
and the stable subspace, respectively) are non trivial. One should actually add an 
analogous condition to the statement of Theorem 16.61 However, to simplify the no- 
tations we can use the convention that [F(0)/((0)] = 0/0 := and that if a subspace 
is trivial then any manifold parameterized by this subspace is just the origin 0. With 
these conventions the statement of Theorem 16.61 is formally correct. We will adopt 
the same convention in the statements of Theorem 16.81 

(5) We now want find an extension of the Uniformly Stable Manifold Theorem. First, 
we recall that Hypothesis H] ensures that there is a manifold of equilibria £ for (16.181) 
which is transversal to the hypersurface {V : ((V) = 0}. One can show that the 
points of S are also equilibria for system (16.161) restricted on the manifold of the slow 
dynamics. 

By relying on the analogy with the toy model, we would like to define a manifold 
A4 US such that any solution lying on Ai us decays exponentially fast to some point in 
S. Also, any orbit lying on A4 US should be written as the sum of a fast dynamic and 
of some stable component coming from the slow dynamics. However, since we are in 
the non linear case we expect some interactions between the different components. 
The precise result is the following. 

Theorem 6.8. Let Hypotheses [IJ . . hold and denote by S the same manifold of 
equilibria as in Hypothesis H Then there exist a constant c > and a manifold 
A4 US which is defined in a small enough neighborhood ofO and satisfies the following 
conditions. 

(a) IfVit) is a solution of (l6~T6D lying on M us and ((V) ^ at t = then ((V) ^ 
for every t > 0. 
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(b) Ai us is locally invariant for (16.161) . Namely, if V e Ai us , then the solution of 
the Cauchy problem 

dt c(vy 1 j 

V(0) = V 

satisfies V(t) e M. us for \t\ small enough. 

(c) IfV(t) lies on M. us then there exists Voo G £ such that 

lim |V(t) - V^o|e c * = 0. 

t— >+oo 

(d) If V(t) lies on M us then 

V(t) = V sl (t) + V f (t) + V p (t), 
where V a \ lies on the manifold of the slow dynamics and Vf satisfies 

lim |^(r)|e CT = 0. 

r— >+oo 

Here the variable r is defined by solving the Cauchy problem (16.191) . Finally, 
V p is an interaction term which is small with respect to the previous two (in the 
sense specified in the statement of [U Theorem 4.2] ). 

The proof of Theorem 16.81 relies on the following lemma of standard flavor. 

Lemma 6.9. Consider the ODE 

^ = F(V) (6.20) 

GST 

and assume that Hypothesis U\ is satisfied and that the center and the stable space 
about the equilibrium point are both non-trivial . Let Sq be a submanifold of a 
center manifold of (I6.20p about the equilibrium point 0. Also, assume that So is 
locally invariant for (I6.20p . Then there exists a constant c > and a manifold M. 
which is defined in a small enough neighborhood of 0, is locally invariant for (16.201) 
and satisfies the following properties. IfV(r) lies on Ai then 

V(r) = V (T) + V f (T) + V p (r), 

where Vq lies on So and Vf satisfies 

lim \V f (r)\e~ = 0. 

T— > + 0O 

Moreover, V p is an interaction term which is small with respect to the previous two 
(in the sense specified in the statement of jH Proposition 3.1}). 



VISCOUS APPROXIMATION VIA ODE ANALYSIS 



31 



Note that, in particular, the statement of Lemma 16.91 ensures that any solution 
lying on M. becomes exponentially close to an orbit on Sq and hence M. can be 
regarded as a slaving manifold for Sq. The proof of Lemma I6T91 is a bit technical, but 
it can be obtained by relying on a fixed point argument. 

Toward a proof of Theorem 16.81 one first defines «So as follows. Consider the 
restriction of (I6.16P to the manifold of the slow dynamics: as pointed out before, 
this system is non singular and £ is a manifold of equilibria. Hence, one can directly 
apply the Uniformly Stable Manifold Theorem 15.21 and define Sq as the uniformly 
stable manifold relative to £ of system (16. 16[) restricted to the manifold of the slow 
dynamics. By applying Lemma 16.91 and by relying on some more work one can then 
conclude the proof of Theorem 16.81 
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